Production Mechanism for Quark Gluon Plasma in Heavy Ion Collisions 



Ambar Jain and V. Ravishankar 
Department of Physics, IIT Kanpur, Kanpur-208 016, INDIA 



PACS: 12.38.Mh, 25.75.-q,24.85.+p 

A general scheme is proposed here to describe the produc- 
tion of semi hard and soft quarks and gluons that form the 
bulk of the plasma in ultra relativistic heavy ion collisions. 
We show how to obtain the production rates in the extended 
phase space - including the colour part - as a function of time 
in a consistent manner, and without having to make ad-hoc 
assumptions. All the required features - the back reaction on 
QCD vacuum, the non-Markovian nature of the production, 
and the quasi particle nature of the partons are naturally in- 
corporated. We illustrate the results with a realistic albeit 
toy model and also show how physically tenable source terms 
may be obtained. 

There is a general consensus that the Quark Gluon 
Plasma (QGP) phase of the Hadronic matter has already 
been produced in ultra relativistic heavy ion collisions 
(URHIC) at CERN. Future measurements at RHIC [1] 
and LHC [2] are expected to provide more information 
on QGP. Be it as may, it is still premature to claim that 
we have a full grasp of the space-time history of QGP. 
Indeed, a considerable effort has gone in addressing the 
issue of evolution of QGP; starting with the pioneering 
work of Bjorken [4] which was itself based on a highly pre- 
scient work of Landau [5] , a number of results have been 
reported [6-8] . All these approaches are based on the idea 
that the bulk of the plasma forms a classical gas, and 
evolves via an appropriate transport equation. Robust 
that this approach is [9,10], it has to be admitted that 
our understanding of the production process is not sat- 
isfactory. Recall that the transport equations (set up by 
incorporating The Yang-Mills (YM) dynamics) require 
source terms that need to be specified in the extended 
phase space (including the compact part corresponding 
to the colour degree) as a function of time; their deriva- 
tion consequently requires care: the basic uncertainty 
principle should not be violated, and the multiple over- 
lapping time scales that pervade the phenomenon should 
not be ignored. They should correctly reflect the dynam- 
ical nature of the vacuum [11], the non-Markovian nature 
of the production rates [12], and the quasi-particle nature 
of the excitations [13]. It is needless to say that a source 
term lacking in one or more of the above of the require- 
ments would predict incorrect and unreliable features for 
the plasma. 

Lattice computations [14] suggest that even at twice 
the transition temperature Tq, the quarks and gluons do 
not behave as ideal gases. We thus need effective field ap- 
proaches based on non-perturbative QCD. The transport 
equation approach mentioned above provides a means of 
introducing collective degrees of freedom appropriate for 



describing the non-perturbative dynamics, via the notion 
of a background chromoelectric field (CEF) which acts as 
a source for the production of the bulk of QGP. The accel- 
eration of the partons in this field takes into account the 
long range interactions amongst them. The short range 
part of the interaction can then be suitably modelled by 
an appropriate collision term. That these aspects can- 
not be captured by the pQCD description [19] has been 
emphasized recently [3]. We thus pay attention to the 
source terms obtained from the CEF degrees of freedom. 

The pioneering work in Ref. [11,12,15] introduces the 
non-Markovian nature into the production process, by 
employing the Schwinger mechanism [17] for particle pro- 
duction. This mechanism is valid only for static elec- 
tric field configurations. Production of particles from a 
space time dependent field would be described perturba- 
tively by the Schwinger Dyson expansion [18]; it is well 
known from transport studies [6,7,16] that the CEF in 
URHIC is necessarily and strongly space-time dependent 
by virtue of energy momentum conservations. See [16] 
for a graphic comparison of the relative contributions of 
the Born and the Schwinger terms. A remedy would be 
to correctly invoke an appropriate perturbative mecha- 
nism [16] (over the nonperturbative vacuum with CEF), 
but unfortunately this study ignores the non-Markovian 
aspect and also makes unjustified assumptions in order 
to extract the time dependence of the source term. Both 
the above attempts completely ignore the essential dy- 
namical nature of the colour charge, which is central to 
the non-abelian plasma that we are interested in. There 
has been some attempt [7] to enlarge the phase space 
to include the colour part; this attempt again involves 
the Schwinger mechanism apart from a lack of rigour in 
the derivation of the colour dependence. Finally, none of 
the above (and other) approaches acknowledges the quasi 
particle nature of the excitations. In short, the source 
terms obtained so far have utmost limited utility, and a 
simple framework containing all the required features is 
still lacking. 

In this paper, we outline a general scheme for deriv- 
ing such a source term, based on simple but general 
considerations, (i) The non-perturbative aspects of the 
mechanism are modelled after the colour flux tube model 
[20] which is known to provide a natural setting for dis- 
cussing quark confinement [21]. (ii) The rate is extracted 
not from an S-matrix computation, but by studying the 
time dependent evolution of the state in the Fock space. 
As a natural byproduct, the quasi particle and the non- 
Markovian nature emerge automatically. 

We now describe the flux tube model in brief. Soon 
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after the two nuclei collide, and start receding from each 
other, colour strings are formed between them. These 
strings merge to form the so called colour rope, a pic- 
turesque way of stating that a CEF is produced between 
the nuclei. The consequent production process tanta- 
mounts to the instability of the QCD vacuum in the 
presence of a classical CEF which is, in general, space 
time dependent. Of course, the dependence itself is to be 
obtained self consistently from a transport equation [6,7]; 
here we merely content ourselves by indicating how to in- 
fer the source term given an explicit history of CEF. We 
study both the qq and the gluon productions in the pres- 
ence of a classical YM field below. We discuss The gluon 
case mainly since its production process is quintessen- 
tially non-abelian with no counterpart in QED. 
Formalism: Consider the gluon case first. The gauge 
potentials A® are expanded as a sum of the classical val- 
ues and their fluctuations 0° :< >= C°. Expand- 
ing the YM Lagrangian, the terms that are responsible 
for the gluon production are identified to be 

L 2g = -§/ o6c [(0„c» - d v ciw h V c 

- ^)(C"V C + r h C vc )] + 0(g 2 ) (1) 

where we have kept terms quadratic in the fluctuations. 
The production rate is to be determined from the above 
Lagrangian, by a suitable choice of the background field 

Unlike in the Maxwell case, the YM configuration is 
not completely specified by the field tensor [24] . It is well 
known [25] that the consequent vacuum polarization also 
differs for each of the above (physically distinct) configu- 
rations. It is thus pertinent to examine which are the con- 
figurations that would make the vacuum unstable against 
particle production. We have checked that a la the re- 
sults for constant CEF [25] only the 'abelian configura- 
tion' of the CEF leads to the quark and gluon productions 
- atleast in the leading order. Keeping in mind the model 
and invoking (local) gauge covariance of CEF, we take the 
form of the gauge potential to be = 5^$ Ci(t, r)5 a ,i 
where the summation is restricted to the diagonal gener- 
ators of the gauge group. For all such configurations, the 
magnetic field vanishes. Denote as usual the initial vac- 
uum state, at t = 0, by |0 >. Let |<5 > (t) = U(t, 0)\0 > 
be the state at any later time t: U (t, 0) is the evolution 
operator determined by the standard Schwinger-Dyson 
expansion that is obtained from Eqn.l. Adapting the 
analysis in [18], the amplitude for the gluon production 
as a function of time is obtained by projecting \*S> > (t) 
to the two gluon state \gg >= \pi,p 2 ',si,s 2 ;ci,c 2 > la- 
belled by momentum, spin and colour quantum numbers 
respectively. In the leading order which we consider for 
simplicity [26], the amplitude is given by 

< 99 \T(W >= • 



C°' a (E 1 +E 2 ;p 1 +p 2 ;t) (2) 
where T(t) = U(t,0) - 1. Further, 

C a - J dt 1 e~ i ^ +E ^ J d 3 r ea ;p(i(pi+|? 2 )-r)C ' a (ti,f) 

is the incomplete Fourier transform of the gauge field and 
Ei are the energies carried by the gluons. The lower limit 
of integration emphasizes that the CEF gets created at 
t = . The corresponding expression for the qq produc- 
tion is given by 

< qq\T(t)\0 >= ^-^=CST^ C2 ul( Pl )v. S2 (-p 2 ), (3) 

with \qq >= \pi,p 2 ; Si, s 2 ; c\, c 2 >. T a are the genera- 
tors of the gauge group in the fundamental representa- 
tion, while u, v are the usual Dirac spinors. For t finite, 
the gluons and the quarks that are in the final state have 
non-overlapping projection with the unphysical degrees 
of freedom. We distinguish the physical pairs by pro- 
jecting the states further onto the physical subspace, by 
imposing the on the mass shell and the transversality re- 
strictions as subsidiary conditions (analogous to the sep- 
aration of the transverse part of the electromagnetic field 
from the longitudinal part in the study of radiation by 
charged particles). 

The probability that a pair is produced any time dur- 
ing the interval (0,t) is given by | < <j>\T(t)\0 > | 2 , 4> 
standing for either gg or qq. The production rate at 
any time t is thus given by its derivative at that in- 
stant; clearly, the rate takes both positive and negative 
values, from which one concludes that the particles ac- 
quire a (time dependent) finite life time. Equivalently, 
vacuum has a dynamical role - both as a source and a 
sink. Further, the rate at any time depends not merely 
on the field configuration at that instant - it has a highly 
non-Markovian character, being determined by the entire 
history. In short, quantum interference effects dominate 
the production mechanism, and this cannot be accessed 
by the approaches based on Schwinger mechanism and 
pQCD. 

A model example: We now illustrate the main re- 
sults with an example drawn from its behaviour expected 
in real situations [4]. Taking the gauge group to be 
SU(3), let the CEF have the form £? = 5^ z £ (8 a . 3 + 
5 a ,&)exp{{\z\ - t)/t }8(t)9(t 2 - z 2 ) with the initial field 
strength given by This choice captures most of the 
features expected of a boost invariant description [4] , al- 
though it is by itself not a boost invariant configuration. 
Notice that the field is characterised by a single time 
scale to, in the units of which all the observables will be 
expressed. Indeed, to gives the decay time of the field 
or equivalently, the overall 'production time' of the par- 
tons, while g is related to the string tension. The initial 
energy density U, to and g are not entirely independent 
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of each other: g U t\ <~ 1. We shall discuss the gluonic 
case in detail. Consider the spin summed rate in the 
momentum space (see below for the colour part descrip- 
tion). Wc display the two particle rates because eventu- 
ally they are the ones needed to study signals involving 
correlations such as dilepton production. Since the rate 
is defined in a six dimensional space, two dimensional 
sections will be plotted as filled contour diagrams after 
scaling them with respect to appropriate powers of to, Eq 
and g 2 . As an example, we fix the total centre of mass 
energy y/s — 0.6 and the transverse momentum pt = 0.1. 
The rates are then exhibited as contour plots in the plane 
defined by time and the longitudinal momentum pi,\ of 
one of the particles. The longitudinal momentum p^2 of 
the other particle then has two solutions, one with the 
same sign as pli, and the other having the opposite sign. 
We designate the two channels by (+,+) and (+,-) re- 
spectively, and show in figs. 1-4 the rates in both the 
channels. Observe the following features: The rate is a 
highly fluctuating function of all the variables involved, 
changing signs rapidly, both with the momenta as well 
as time, as shown in Figs 1-4. The rate oscillates over a 
much larger range for smaller values of pr at a given y/s. 
For example, it oscillates roughly between (—500, 500) at 
y/s = 0.6, pt = 0.1 (Figs. 1, 2). In contrast, the range 
is restricted to <~ (—100, 100) when pt is increased to 0.2 
(Figs. 3, 4). Notice also that for a given y/s and pr, the 
fluctuations are more dominant for very small and very 
large values of longitudinal momenta. We have checked 
that the rate and its fluctuations are larger at small val- 
ues of y/s, dropping rapidly as y/s is increased. Though 
not shown here, the number of oscillations are also seen 
to increase as we increase y/s. It may also be seen from 
the figures that the (+,+) channel is relatively placid 
in comparison with the (+,-) channel. However, the ac- 
tivity at intermediate values of longitudinal momenta is 
more pronounced in the (+,-) channel for a given value of 
y/s,px- In short, these figures exhibit in a vivid manner 
the production and the the absorption of the pairs, and 
the correlation between them, in the momentum space. 
The correlations are of vital importance in getting the 
screening lengths as well as in studying the hadroniza- 
tion. Coming to the quark production, suffice it to men- 
tion that they are equally rich in features differing only in 
the topology of the distributions. We do not show them 
here. 

The analysis is incomplete in that the rates are still not 
described in the phase space. This has been a vexed issue, 
with solutions attempted by imposing additional condi- 
tions rather arbitrarily [4,27,6]. We resolve this issue by 
projecting the amplitudes to a basis of coherent states in 
the extended phase space, a method that allows as accu- 
rate a classical description as possible for these quantum 
processes. As an illustration, consider the colour part 
which is unique and central to the plasma under dis- 
cussion. The structure of Eqns.(2,3) guarantees that it 



is independent of momentum and spatial variables. Let 
the gauge group be SU (2) for simplicity (it has a simple 
topology). Denoting the two gluon colour state in the 
coherent basis by \h\,fi2 > [28], the colour part of the 
source term may be written as Y, 9 (ni,fi2) ~ \n\ x tt-2 | 2 
This expression has been obtained by averaging over the 
direction of the CEF in the group space. Contrast S 9 
with its counterpart for qq production which has the 
structure Y, q (hi,h 2 ) ~ \{l — hi ■ n 2 ). While the glu- 
ons are predominantly produced in directions normal to 
each other in the colour space, the quark pair prefers 
to be produced in antiparallel directions, a feature that 
holds for SU(3) as well, although in a more involved fash- 
ion due to its more complicated topology. Coupled with 
the fact that the colour charges are indeed vectors in the 
group space and that they are dynamical variables obey- 
ing the Wong equation [29], we conclude that the dyan- 
mics of QGP would differ significantly in the bosonic and 
fcrmionic sectors in a manner that cannot be captured by 
simply multiplying the Maxwell results by the number of 
degrees of freedom. Finally, we have performed a pre- 
liminary computation and found that the phase space 
dependences of quarks and gluons in position and mo- 
menta bear no resemblence to each other; Those details 
will be presented in a separate communication. 
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